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Abstract 

O ' We derive the random-phase approximation for spin excitations in general 

multi-band Hubbard models, starting from a collinear ferromagnetic Hartree- 
^ ' Fock ground state. The results are compared with those of a recently intro- 

o 

f^ ' nets. As we exemplify for Hubbard models with one and two bands, the two 



duced variational many-body approach to spin-waves in itinerant ferromag- 



approaches lead to qualitatively different results. The discrepancies can be 
traced back to the fact that the Hartree-Fock theory fails to describe properly 



n^ • the local moments which naturally arise in a correlated-electron theory. 
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I. INTRODUCTION 

Metallic ferromagnetism is one of the oldest problems in solid-state physics. Despite all 
progress in solid-state theory over the last centuryETH, we are still far away from a definite 
theory for transition metals like iron, cobalt and nickel. It is generally accepted that the 
magnetic order in metals is a consequence of the interplay between the electrons' kinetic 
energy and their mutual Coulomb interactions. Therefore, it does not come as a surprise 
that a complete theory for this many-particle effect is difficult to develop. 

Nevertheless, it is commonly believed that at least the ground-state properties of itin- 
erant magnets can be described by (effective) one-particle theories such as Hartree-Fock 
theory or (spin-)density-functional theory (SDFT). However, serious doubts are in order, as 
a closer look on available experimental data reveals. For example, SDFT fails to reproduce 
even qualitatively the electronic band-structure of nickel, one of the simplest itinerant fer- 



romagnetsu; for a summary, see Ref. |jTO|]. For iron, the agreement is better but cannot be 
considered as satisfactory either. In any case, the same physical mechanism is acting in both 
systems, and a convincing theory should cover all elements of the iron group equally well. 

Single-particle theories are not very satisfactory also from a purely theoretical point 
of view: It is well known that they significantly overestimate the stability of ferromag- 
netismtl. For example, Hartree-Fock theory for one-band Hubbard models gives a ferro- 
magnetic ground state above some moderate interaction strength. In contrast, many-body 
theories show that ferromagnetism in one-band Hubbard models occurs, if at all, only for very 
large Coulomb interactions or very peculiar choices of the one-particle density of stateslij'li3. 
This indicates that the Stoner mechanism is less prominent in itinerant ferromagnets than 
implied by single-particle theories. 

As already argued by van VleckEl and confirmed recentlycjilj, an essential ingredient for 
metallic ferromagnetism is the local (Hund's-rule) exchange coupling of electrons, a generic 
feature of electrons in degenerate bands. In contrast to the one-band case, ferromagnetism 
occurs for moderate Coulomb-interactions in a generic two-band model provided that local 



Hund's-rule terms are included. Therefore, even a qualitative understanding of the 3d- 
elements requires a careful consideration of their atomic open-shell structures. 

In our recent worktJO we have developed a variational approach which allows us to exam- 
ine multi-band Hubbard models with an arbitrary number of correlated orbitals per lattice 
site. For a realistic description of the Sd-elements a minimal model includes the 3d, As, and 
Ap bands. When applied to nickel, our theory provides results which agree much better with 
experiments than all previous SDFT-calculationli2il3. In particular, our approach reproduces 
the experimental Fermi-surface topology of nickel and removes all other qualitative SDFT 
deficiencies. Since our variational many-body approach provides a reasonable description of 
the ground state, it may serve as an appropriate starting point for the study of excitations 
in itinerant ferromagnets. 

As observed in inelastic neutron scattering experiment scJ'E3, the spectrum of itinerant 
ferromagnets displays low-lying spin-wave excitations. These excitations are generally con- 
sidered to be the driving force for the transition from the ferromagnetic phase into the 
paramagnetic phase at the Curie temperature. Therefore, their correct description is also 
very important for our understanding of finite-temperature properties of itinerant ferro- 
magnets. Unfortunately, almost all established theories which are used for the analysis of 
spin-waves in transition metals are based on effective single-particle theories which do not 
necessarily provide a satisfactory description of the ground state. 

The textbook approach to the problem of itinerant ferromagnetism is a Hartree-Fock 
theory for the ground state which is then combined with a random-phase approximation 
(RPA) for the description of spin-excitationg^u. For example, in Ref. |^ the RPA-method 



has been applied to iron and nickel starting from a SDFT calculation for the ground state. 
Later, a similar method has been used to analyze data from electron energy-loss spectroscopy 
for the same materials£3. Another common approach starts from energy calculations for 
static spin deformations of the ferromagnetic ground state within the SDFT. The spin-wave 
energies are identified as the energy difference between these "frozen-magnon" states and 
the ground states. Despite their conceptual shortcomings both SDFT-based approaches are 



able to reproduce the experimental data surprisingly well. 

Very recentlycJ one of us has proposed a spin-wave theory which starts from the varia- 
tional ground state of multi-band Hubbard models as introduced in Ref. [0. This theory 
allows us to determine the spin-wave spectrum for systems with an arbitrary number of 
correlated orbitals per lattice site. The main advantage of this new approach lies in its foun- 
dation on a true many-particle description of the ground state. In this way, we are able to 
provide a consistent picture of the ferromagnetic Fermi liquid and its spin-wave excitations. 

In order to obtain a better understanding of our own but likewise of the SDFT-based 
methods it is crucial to apply them for a generic but simple model system. To our knowledge 
such an analysis has not been performed yet. Even a comparison between the two well- 
established SDFT methods is still lacking. In this work we analyze the Hartree-Fock/RPA 
theory and compare it with our variational many-body approach. Although the RPA is 
generally considered as the "standard" method for the description of magnetic excitations 
in itinerant ferromagnets [Q there exists no systematic derivation for multi-band Hubbard 
models. In Sect. fT|, we derive the general RPA equations from a decoupling scheme and a 
standard diagrammatic viewpoint. 

In Sect. |T| we provide numerical results for a generic two-band model and show that both 
methods lead to qualitatively different results. For generic values of the magnetization, the 
Hartree-Fock/RPA treatment incorrectly predicts an instability of the coUinear ferromag- 
netic phase against some spiral ordering. In our variational approach the spin-waves always 
have positive energy and their stiffness decreases with increasing correlations. The case of a 
one-band model with fully polarized bands allows us to pinpoint the origin of these discrep- 
ancies: The Hartree-Fock theory allows for ferromagnetism with small local moments and 
large longitudinal spin-fluctuations due to Stoner excitations. These are actually absent in 
a realistic correlated-electron description where charge fluctuations are strongly suppressed 
in the ferromagnetic and paramagnetic phases. 

Conclusions in Sect. |IV| close our presentation. 



II. RANDOM PHASE APPROXIMATION FOR MULTI-BAND SYSTEMS 

In this section we derive the random phase approximation (RPA) for ferromagnetism 
in multi-band Hubbard models. Starting from a LDA-calculation, RPA results for iron 
and nickel have been reported in Ref. [^. However, in these calculations the relevant 
Coulomb matrix elements were significantly simplified. In this work we present a more 
general treatment. 

A. General remarks 
We will address the following general class of multi-band Hubbard modelsll3'E3 

^ = E E ^M^sc,,s' + E ^;- = ^1 + ^t . (1) 

i,j s,s' i 

Here, cf^ creates an electron with combined spin-orbit index s = {b, a) at the lattice site Ri 
of a solid {b = 1, . . . ,N] o"=t,|; A^ = 5 for 3d electrons). The atomic Hamiltonian 

Sl,S2,S3,S4 

is assumed to have site-independent interaction parameters ^^i'**2;s3,s4_ 

In order to determine the spin-wave properties, we need to consider the imaginary part 
XriQ) E) of the transverse susceptibihtyiS, which is given by the retarded two-particle Green 
function 

GUq,E) = ^{{Sl,S-.))E (3a) 



L.,0 



y""dte^^*(vl/o|[^+(t),5^::(0)]|vl>o) . (3b) 

Here, we introduced the g-dependent spin-flip operators 

/ l,b 

s,=isp^ = Y.^-''%Hfm)^ (4b) 



Lb 



in the Heisenberg picture, where the sum runs over all L lattice sites Ri and orbitals h. 

Magnetic excitations are found as poles of the Green function Gret(^, E)^ or, equivalently, 
as peaks in 

XT{q.E) = luiGU^,E) (5) 

at energies E > {). The Lehmann representation of (^, 



(^n I Sz I M/o) p I (vj/, I 5+ I v&o) 



E - {En -Eo)+tr] E + (E„ - Eq) + tv 



(6) 



shows that the weights of the poles in Cretin, E) at the energies En — Eq > are given by 

Wn =1 i^n I 5-1 I vPo) |2 . (7) 

In a ferromagnetic system the state 

i^;-) = ^,=^i^o) (8) 

is also a ground state of if for g = because the operator SZ_^ just flips a spin in the 
spin multiplet of the ground state |\l'o)- Therefore, we can conclude that Gt(0, E) has one 
isolated pole for E — Eq = 0. This is not surprising since the spin- wave is the gapless 
Goldstone mode in the symmetry-broken ferromagnetic phase. For finite values of q, it is an 
experimental fact that there are also pronounced peaks in Xriq, E) at the spin- wave energies 

E = E{q) ^ Dq^ . (9) 

The constant D is usually denoted as the "spin-wave stiffness" . 

B. Hartree-Fock treatment 

The one-particle Hamilton Hi in (0) may be written in momentum space as 

k;b,b',(T 



fc;«,o',o- 
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Here, we introduced the operators 

N 

6=1 

AT 



k:a,( 



k:a,: 



(11a) 
(lib) 



6=1 



and the respective energy matrix-elements 



Eaa' {k, (^) = Y1 ^b,b' {k)u*^/k, a)M„, ,,, {k, a) , (12) 

6,6' 

where the elements Ua,b{k, cr) of the unitary matrix u{k, a) will be identified later. The 
transformation (|Tl]) leads to the following expression for the atomic Hamiltonian, 

(13) 
where 

^Iflf iPl,P2, P3,Pi) = Yl ^a,a'''^''^'«^i,6i (PU (^)Up^,b2 (P2, Cr')Up.^^bs (PS, f^')W/34,64 (^4, Ct) , 

61,62,63,64 

(14a) 



7-7-61,62,63,64 _ 7/(61 o-),(620-');(630-'),(64ct) 



(14b) 



We choose the matrix u{k, a) in order to fulfill the corresponding Hartree-Fock equations 



Ef{k, a)5pp, = Ep^ik, a) + 2 Yl f/^^^"""'^ [<,64^;,62(^, ^K'^b^ii cr) 

61,62,63,64,0-' 

-S:'T,^;^,.yf,^,^{k,a)Uf,,^,^{k,a) 



(15) 



with 



'^b,b' = 2^1 ^?^h(P^ 0^)^/3,6' (P> 0^) 

P,/3 



n: 



f3a 



^ 



HF 



^p;l3,a(^p;f3,a 



^ 



HF 



(16a) 
(16b) 



Our Hartree-Fock theory is restricted to cases where the translational invariance is conserved 
but we still allow for a coUinear ferromagnetic ground state. The Hartree-Fock Hamiltonian 
thus reads 

^"" = E ^r(p»«;AA;/3,. • (17) 

In the special case where, (i), all orbitals h belong to different representations of the point- 
symmetry group of the lattice, and, (ii), the Hartree-Fock ground state I^E'^^) is invariant 
under the respective symmetry operations, Eq. ( [16a| ) becomes 

^M' = h^nt ■ (18) 



Under these conditions, the results of Ref. ^^ are recovered; see below. 



C. RPA from the equation-of-motion technique 

We will derive the Green function (^) in the random-phase approximation using the 
standard equation-of-motion techniqueE^. First, we rewrite Gret{<i-,E) as 

GUlE) = Y,G'^et\lE), (19a) 

b,b' 

G\t'\l E)^\ E((^t.1..T%..i' %'M,,'a))- ■ (19b) 



L 

k,k' 



The equation of motion for Gj-et^' i'l-iE) cannot be decoupled self- consistently. Therefore, 
we will consider the more general Green function 

which allows us to express ( |19b| ) as 

Gi^yiiE) = 1 5^ E<,^i(^+^"''T)««...(^'i)E^ltS(^) • (21) 

ai,Q2 k k' 

The equation of motion for G^_^^^y{E) has the form 



EG'^^;^l,,{E)=A + B + C, (22a) 

^ = ^^^t,^auA,^2,l ' ^L'/k'-,r,b',M^ ' (22b) 

^ = <a^J.,«.T^M.i' ^-] ' 4;.a%.^;..T))- • (22d) 

In RPA we assume that the exact ground state |\l/o) of our Haniiltonian (|l]) may be replaced 

by the (spin-polarized) Hartree-Fock ground state 1^1/^^) of (p^ . Then, the expectation 
value ( p2b|) becomes 

A = hrk'<Ak^ Du^Mf^+q, T) (^^;;- ^r') • (23) 

The evaluation of (|22c|) leads to 

B = Y1 ^--/^(^' i)^Sf.'(^) - E ^/^.- (^ + ^"'' t)5|B'(^) • (24) 



Finally, the commutator in ( p2d| ) generates two-particle operators, which we decouple ac- 
cording to the rule 

((0+0+03^04^; . . .))e ^ (a^|a3-f)vtHF((a+a4^; . . .))e - {a+^a^^)^HF{{a+^d^^; . . .))e ■ (25) 

This approximation appears to be natural because we replaced |\I'o) by the one-particle 
product state |\I^[^^). Note, however, the replacement (|25|) does not become exact even if 
we work with 1^^^) instead of l^&o)- 

After the application of the decoupling scheme ([25|), Eq. (p2d|) may be written as 



C = Ci + C2 + C3 (26a) 

2 



Ci = ^Y.Il [<V,t'' iP^ k, k, p) - <.V-f {p, k, p, k )] Gf;% (E) , (26b) 

/3,/3' p;a 

(26c) 
^3 = I E E (^r' - ^£',) ^;A? (p + <l^ k, p, k + g )Gg;|, (E) . (26d) 



When we use the exphcit expression ([T^), the sum over p can be carried out for Ci and C2 
because the Green functions do not depend on p. 

c, = 2 Y. u''^'^'"'' E ^i' [(^i^i,^4 + Tl,A) <,i>.ik, i)u,^,M i) 

SiTUu:,,,,{k, l)u,^,^{k, i)] Gf'5f,,(i?) , (27a) 

Tl^bA'^P^i^ + ^^ T)Mai,fe4(^ + ^> T) 



bi,b2,b'^,b4,cr,a' 



C2 = 2 Yl ^:.'"'"'^E^i^aS'(^) 

bi,b2,b3,b4,,a;a' /3 

%tU + ^i'^U ) «A^i(^ + <i^ T)«a„.,(^ + q, T) 



(27b) 



^3 = -2 E f/S'^^'^^'^^(n^;;-nf)<,^(fc;i)«,^,^(fc>g,T) 



biMMM 



Yl ^ftfel(^+^'T)W/3^63(^'i)<^S';'^>'(^) • 



(27c) 



p;/3,/3' 



With the help of the Hartree-Fock equations ( [T5| ) and relation (|2TD we may cast the Green 

function ( [20| ) into the form 

-1 



yk,k' ,q 



-.HF/ 



^^^^..(i?) = i? - < (fc, i) - i?„",^ (k + q, T) 



%,a:'<,.'(^, i>a,Ak + q^ T) (^^;i- ^?' 



(28) 



-2 5^ G!:^^''''(g,E)f/J|' 
This equation, together with (pT]), leads to 



^2,^3,^4 /^QlT ^C»2i 



\^;.-- ^? j <,fe2(^' i)^.„6,(/^ + g, T) 



G'il^''''{q, E) = [rret(g, i5)^,t] + E [rret(g, i?)££] uS^G'iiyiq, e) 

b3,b4,b5,b(i 



(29) 



where 



r,M,E)lili^^ Y 



<i,6i(^ + ^,T)Ma2,fe2(^, i)Mai,63(^+^'T)<2,64(^'i). „i| „2| 



k,ai,a2 



Tjb-ibi _ _cyTTb5,bi,b(,,b-j, 



E- EHF(^J)_^HF(fc + ^-.J) +,^ 



n^ ^— ri- 

k+q k 



(30a) 
(30b) 



Here, we added the infinitesimal increment rj = 0^ to ensure the properties of a retarded 
Green functioned. When we consider Tj-et{q,E) and U in ( PPJ ) as matrices with respect to 
the indices 2) = (6162)5 the solution of Eq. (^) is given by 



10 



Gf/(g, E) = J2\{^- rrct(g, E) ■ u 



-1 



3' 



rret(g,i5)^,, . (31) 



JS' 



This result, together with Eq. ( |19a|) , gives us the Green function Gret{(l,E), whose poles at 
E^^^{q) define the RPA spin- wave dispersion. We will further analyze it in section |IT1[ 

D. RPA from the diagrammatic approach 

The result (^T]) of the RPA decoupling scheme may also be derived from the standard 
diagrammatic approach. Here, the matrix ( |30a|) is nothing but the retarded part of the 
respective causal transversal susceptibility F, 

^lllM E) = ^Y.j dE'Gfl,{q + k,E + E')Gfl^,{k, E') , (32) 

k 

where the one-particle lines are evaluated within the Hartree-Fock approximation; see be- 
low. Using this expression and the vertex for spin-fiip excitations t^f^^ from ( |30b|) , we can 
calculate the causal two-particle Green function G^' (g, E) as the usual RPA sum of the 
"bare bubbles" in (BSI), 



G^^'\q, E) = Vf,M E)+Y1 ri'(?"*, E)Ug,rZ{q, E) + . . . (33a) 

XI', D" 

= Tf,,,% E)+J2 ri'(g, E)US:G^"'''{q, E) . (33b) 

Xi',D" 

This equation has the same solution as (|3T|) if we replace Fret and Grct by F and Gt, 
respectively. The retarded Green function Gj-et can then be re-derived from Gt with the 
help of the standard relationsH^I 

ImGret(g,^) = lmGT{q,E) {e{E) - &{-E)) , (34a) 

Re G,et{q, E) = Re Griq, E) . (34b) 



In order to show the equivalence of both approaches, it remains to calculate F in 
explicitly. To this end we first determine the Hartree-Fock Green function Gf,^^{k,E). Its 
diagrammatic evaluation leads to 

11 



G«^,(fc, E) = Gl^^ik, E) + J2 G%,AK E)T}^l,^„G'tlSK E) 



(35) 



b'h" 



where G^y ^{k, E) is the causal Green function corresponding to the one-particle operator Hi 
in®, 



Glik.E) = E5t,b' - et,Ak) . 

J bb' 



Moreover, we introduced the usual Hartree-Fock self-energy 



(36) 



Y'HF 9 \ ^ I TTb,bi,b2,b' rpcr' i:ct' Tj-bfii.b' ,b2rpa 

feib2,o"' 



(37) 



whose elements are calculated self-consistently, 



Tb,b' = E (^ 



'T 



^k;(bu)%{b'a) 



^r 



(38) 



compare ([T^). Using a matrix notation in ( P^D we find 



Gf (fc,E) " = G°(^,E) 



^HF 



E-E,{k) 



(39) 



where the matrix E(j{k) is given by 



i5;.(A:)) =e,,,,(fc) + E^,^,^. 



6,6' 



(40) 



From Eqs. (l39D and ( ^OD we see that [(^^^(A;, -E)]"^ is diagonalized by the matrix u{k,a) 
which solves the Hartree-Fock equations ([T5|). The eigenvalues are given by E^^{k, a). The 
diagonalization and inversion of (|39D give the final result 



1 



G^'{k,E)] =6sB'- 

'PP' ^'^ E-Ef^k,a)+ir]p{k,a) 



(41) 



Again, we added a positive or negative increment rjjs^k, a) depending on whether E^^{k, a) 
is larger or smaller than the Fermi energy. This ensures the analytical properties of a causal 



Green function. We insert fETf) into 



and find 
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Tlf^{q,E) = \ Y. <,6,(^ + ?"',TK„..(^,iK,.3(^ + ?"''TK„6,(^,i) 



L _ 



1 — n^ jn-. 



X 



( -, ^±^^ r (42) 

^E - [E^^{k + q, T) - E^lik, i)j + zri 






E - [E^fik + g, T -E^!{K i))) - ^r 

From this equation one easily re-derives (|30a|) with the help of the general relations (p4D. In 
this way, we have shown the complete equivalence of the diagrammatic and the equation-of- 
motion derivation of the multi-band RPA equations. 

III. SPIN- WAVE DISPERSIONS 

A. Variational spin-wave dispersion 

In Ref. [|^ we proposed the following Gutzwiller wave-functionS for a variational exam- 
ination of the Hamiltonian (|I|): 

\^g) = Pg\%). (43) 

Here, |$o) is any normalized single-particle product-state and the Gutzwiller correlator Pq 
allows for a variational adjustment of the occupation of local atomic multiplets. Our choice 
for the correlator Pq ensures that the wave-function (^) yields the exact ground state of H 
both in the uncorrelated and the atomic limit. For all other values of correlation parameters 
the expectation value of H in the wave-function can be determined analytically in the limit 
of large spatial dimensions, d -^ 00. As in our recent work, this limit must be considered 
as yet another approximation since we apply our general analytical results to real three- 
dimensional systems. Note that l/tZ-corrections are expected to be smallli3EZl. 

Once the optimum variational ground state \'^°q ) has been found by minimizing the 
ground-state energy-functional, the following expression for the variational spin-wave dis- 
persion can be evaluated 

13 



E^^\q\ 



^T 



SiHSz 



^opt\ 



^opt 



sp, 



^opt 






(44) 



In Ref. [^ a general analytical expression for E^^^{q) has been found in the limit of large 
spatial dimensions. The numerical evaluation of this result for real materials like iron or 
nickel is involved but feasible. So far, we have applied our method only to a degenerate 
two-band model. The results for this model will be compared with the RPA in the next 
subsection. 



B. Two-band model 



In Ref. [P3| we have discussed the variational spin- wave dispersion for a model with two 
degenerate Cg orbitals (6 = 1, 2) on a simple-cubic lattice. In this system the general atomic 
Hamiltonian (^ becomes 

Hs,t = U^ nb,^nb,i + f^' ^ ^i,<7^2,<7' - J^ ^i,<7^2,<7 (45) 

b cr,a' (J 

+ J 2_^^l,a'^2,-(7'^l,-a^2,a + '^C I Cj^ ^C-^ ^Cg ^Cg | "I" ^2 ^Cg |C]^ ^C-^ ^ 1 . 
u 

In cubic symmetry, the Coulomb and exchange-integrals f/, f/', J and Jq are not independent 
from each other. Instead we have only two free parameters, because the relations J = Jq 
and U-U' = 2J hold. 



In Ref. [0] we have discussed the appearance of ferromagnetic order in this model in 
detail both for the Hartree-Fock and our Gutzwiller theory. As is well known for mean-field 
theories, the stability of a ferromagnetic solution is drastically overestimated in a Hartree- 
Fock treatment. In a Hartree-Fock description only the Stoner parameter I = {U + J)/2 
governs the magnetic behavior. The paramagnet becomes unstable if the Stoner-criterion 
I Nep > 1 is fulfilled, where Ne^ is the density of states at the Fermi level. This means that 
a ferromagnetic transition occurs for any value of the local exchange constant J. This is 
in striking contrast to our variational approach. In the Gutzwiller many-body approach it 
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becomes evident that a sizeable Hund's-rule coupling J is crucial for the formation of ferro- 
magnetic order. Thus, for small values of J, the Hartree-Fock theory leads to qualitatively 
incorrect results. 

The significant differences for the critical values of the Coulomb interaction complicate 
a comparison between the RPA and the variational spin-wave approach. For the same set of 
parameters J and U the underlying ground states are completely different and, therefore, a 
comparison of the spin-wave properties would not make sense. Thus, it appears to be more 
reasonable to consider the results of both methods for the same magnetization per band, 

^ = \ H^S-Q = ^ E^^'^^t) - (^/;h) • (46) 

I I 

The density of electrons per band and spin direction in the paramagnetic phase is given by 

In Fig. 1^ the variational spin-wave dispersion in x-direction is shown for four different 
magnetizations. Here we used the same tight-binding parameters as in our analysis of the 
ground-state properties in |TB|. The average electron density per orbital and spin direction in 
the paramagnet is n = 0.3. We keep the ratio J/U = 0.2 fixed and consider different values 
U/eV = 7.8, 10, 12, 13.6 (bandwidth W = 6.6 eV^) which correspond to a magnetization 
per band oi m = 0.12, 0.20, 0.26, 0.28. The last value belongs to the almost fully polarized 
ferromagnet {m ^ n). As can be seen from Fig. |I], the spin- wave dispersion drastically 
depends on the magnetization. The data can be fitted very well to 

E^^^iiq,, 0, 0)) = Dqlil + ^ql) + 0{ql) , (47) 

in qualitative agreement with experiments on iron-group metals. In the case of a strong 
ferromagnet, we obtain D = 1.4 eVA^ and D = 1.2 eVA^ for m = 0.26 and m = 0.28, 
respectively. This is the right order of magnitude in nickel where D = 0.43 eVA^. The 
inset of Fig. |1| shows that the spin-wave dispersion is almost isotropic which also agrees 
with experimental observations. In our variational many-body approach, the spin-wave 
stiffness decreases as a function of U. This does not come as a surprise because the effective 
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coupling between the sites decreases when the electron hopping between the sites becomes 
less effective. 

In Fig. ^ the results for the spin-wave dispersion in the RPA are shown for three different 
values of U. In contrast to our variational approach, the RPA predicts a negative slope of 
the dispersion for non-saturated ferromagnetism, m < n. Note that for U = 5.2 eV the 
system has just attained the fully polarized state (compare Fig. 7 of Ref. ||l3l). But even in 



this case, a well-defined RPA spin- wave excitation does not yet exist. It requires a further 
increase of U until the fully polarized coUinear Hartree-Fock ground state appears to be 
stable. Since the RPA probes the 'local' stability of a mean-field state it is conceivable that 
another ferromagnetic Hartree-Fock state with some spiral order has a lower energy than the 
collinear state. However, ferromagnetic phases with spiral order are not generic for transition 
metals. We argue that this prediction of exotic ferromagnetic spin structures is an artifact 
of the Hartree-Fock approach. In fact, as we have demonstrated for our two-band modeled, 
Hartree-Fock theory predicts ferromagnetism in regions of the {U, J) phase diagram where 
the paramagnetic phase is actually stable. 

In the case of a stable collinear ferromagnetic Hartree-Fock ground state, the RPA spin- 
wave stiffness increases as a function of U, again in contrast to our variational approach. For 
all generic cases, D^^^{U -^ oo) > D^^^{U -^ oo) holds in the limit of infinite coupling. In 
our case, D^^^{U ^ oo) = 1.4 eVA^, as compared to D'^'^iU ^ oo) = 0.4 eVA^. The vari- 
ational and RPA results differ significantly in the strong-coupling limit, too. Nevertheless, 
there are some intermediate coupling strengths f/o, Jq where the variational and the RPA 
spin-wave stiffness agree. Such accidental agreements may also occur when RPA results are 
compared to experimental data. 

In some special cases where the fully polarized state is a one-particle product state we find 
D^P^(f/, J) < D™(f/, J) for all finite values of f/, J, and D^^^{U ^ oo) ^ L)™(f/ -^ oo) 
for large couplings. This is also the case in the one-band system; see below. 
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C. One-band model 

The study of the one-band model allows us to elucidate the origin of the qualitatively 
different results for the spin-wave dispersion. In the case of the one-band model, 

^=J2Y1 ^i'^^t,aCj;a + ^Y1 ^^T^ii = Yl ^(^) ^ta%a ^ ^ Yl ^*T^*i ' (48) 

i,j cr,(T' * fc * 

the results of Sect. |T^ simplify significantly. The Hartree-Fock energies can be written as 

E^^{k, a) = e{k) + Un_„ (49) 

and the susceptibility ( |30a| ) becomes 

r,et(g, e) = \Y. P'f ~ ""? A ' (50a) 

^ g E-/\- [e{k) - e{k + q)j + ir] 
A = U{n^ - ni) . (50b) 

For the spin- wave energies E^^^{q) the denominator in 
vanishes, i.e., 

i + f/r,et(g,i5;«^^(g)) = o. (52) 

To the leading (second) order in g we can expand (|50a]) in powers of 1/A because 



\E - {s{k) ~ sik + q))\ <^ A (53) 

is fulfilled for our low- lying spin- wave excitations, see (|). This expansion leads to 
^^"^(^l = ^^^^-T^ [E(%^..-T - -h) {<k) - eik + g )) 

k 
k 

This result is to be compared with our variational spin-wave dispersion. However, such a 
comparison is somewhat questionable, since the existence of a stable ferromagnetic phase is 

17 



hard to obtain in our variational approach. Nevertheless, let us assume that a given one- 
band model with a special density of states may have a fully-polarized ferromagnetic ground 
state. Then we obtain the following expression for the variational spin-wave dispersion 

^™^(g ) = ^Y. ""U,! (^(^) - ^(k + q)) . (55) 

k 

This expression differs from E^^^{q) in ( |54D by the second term, i.e., for rij = we find 

5E = E^P^(g)--E™^(g) = -l-LX^n,^^(5(^)-.(fc + g))' . (56) 

k 

The energy difference 6E is always negative and vanishes for U -^ oo, i.e., for A -^ oo. It is 
essentially this contribution which eventually leads to negative values of the RPA spin-wave 
stiffness for finite U in the two-band model. 

From a formal point of view the origin of the differences between both theories is obvious: 
For a fully polarized ground state of the one-band model the variational dispersion is identical 
to the exact first moment of the spectral function XTet{q,E) = lm.Gj-et{q,E). However, 
the RPA predicts the existence of spectral weight around E = A, the so-called "Stoner- 
excitations" . Thus, the spin- wave pole of the RPA must have a lower energy than E'''^^{q), 
since, otherwise, the RPA could not give the correct first moment of xriq^E). For this 
reason we need to address the issue of whether or not the RPA-theory correctly describes 
the high-energy physics of itinerant ferromagnets. 

First, it should be noted that the existence of the energy scale A is an artifact of the 
Hartree-Fock theory which survives in the RPA. From ground-state concepts like the mag- 
netic condensation energy we infer that such an energy scale is an artificial feature of mean- 
field theories, and irrelevant within a many particle description of itinerant ferromagnetsEJ. 
Second, the misleading relevance of this energy scale stems from the overestimation of lon- 
gitudinal fluctuations in mean-field theories. 

In order to see the second point, we consider again the spectral function XTiq,E), (||), 
which provides the energy distribution of the spin excitations. In a spin-wave state the local 
moments do not depend on momentum, i.e, {Si ),j,o = {Si )^o. In particular, this implies 



for the fully polarized one-band model that there are no doubly occupied lattice sites in 
|\l/°.). The temporal development of this state is now crucial: In a strongly correlated 
many-particle system, longitudinal fluctuations are substantially suppressed because charge 
fluctuations are energetically too costly at large U. Therefore, the length of local spins {Si ) 
is basically conserved also as a function of time. Our variational method correctly captures 
this generic behavior: fast longitudinal fluctuations are small and may be neglected at least 
for small values of q. In contrast, longitudinal fluctuations are always present in mean-field 
theory since electrons with opposite spin meet each other with non-zero probability. Recall 
that Hartree-Fock theory predicts a ferromagnetic ground state for much smaller interaction 
strengths than a proper many-body treatment. The presence and importance of considerable 
longitudinal spin fluctuations within Hartree-Fock and RPA then leads to the prediction of 
an instability of the collinear ferromagnetic phase {D^^^ < 0), as seen in Fig. ^ 

These findings also explain the qualitatively different behavior of the spin-wave stiffness 
as a function of U. In the variational approach we find a stable ferromagnetic ground state 
with well-defined local moments. The interaction of these 'spins' is mediated by hopping 
processes. Similar to the well-known large-f/ antiferromagnetic coupling Jaf oc t'^/U in 
the case of the half-filled Hubbard model, we find that the effective spin-coupling Jp in 
ferromagnets is reduced with increasing U because hopping processes are more and more 
suppressed as they induce charge fluctuations. Note that the coupling Jp remains finite 
in the limit U ^ oo because we are dealing with a non-integer band-filling where electron 
transfers are never completely suppressed. In contrast to this generic behavior, the Hartree- 
Fock theory misleadingly predicts a ferromagnetic ground state already for rather small 
interaction strengths where local moments are only partially generated. The RPA dynamics 
of the 'Hartree-Fock spins' is considerably determined by longitudinal fluctuations. With 
increasing interaction strength the local moments are stabilized even within Hartree-Fock 
theory. Consequently, the RPA spin-wave stiffness increases and eventually overshoots the 
value from our variational approach. The basic mechanism for a reduction of Jp in our 
variational description remains completely hidden in the RPA because the reduction of 
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the electron transfer amplitudes due to the electrons' mutual Coulomb interaction is not 
contained in a Hartree-Fock mean-field approach. 

IV. CONCLUSIONS 

In this work we have derived the random-phase approximation (RPA) for the analysis of 
spin-wave properties in general multi-band Hubbard models with a collinear ferromagnetic 
ground state. The equation-of-motion technique and the diagrammatic approach lead to 
identical expressions. 

We have applied our analytical results to model systems with one and two correlated 
orbitals per lattice site. This allows us to compare the RPA results with those of a recently in- 
troduced variational approach to spin-wave excitations of Gutzwiller-correlated many-body 
ground states. The numerical analysis of the two-band model reveals significant qualita- 
tive differences between both methods. In the RPA, the spin-wave stiffness is negative for 
non-saturated ferromagnetism. With increasing correlations the RPA spin-wave energies 
are shifted upwards and the RPA spin-wave stiffness eventually becomes positive in the 
saturated ferromagnet. Our variational correlated-electron approach shows a completely 
different behavior. The spin-wave energies are always positive and the spin-wave stiffness 
decreases with increasing interactions. The same discrepancies occur in the one-band case 
where a fully analytical evaluation of the spin-wave dispersions is possible. 

We conclude that the Hartree-Fock/RPA theory provides an inadequate description for 
spin-wave properties in correlated electron systems. In our view, the main problem is the 
mean-field character of the underlying Hartree-Fock theory. In this approach the ferromag- 
netic ground state lacks well-developed local moments, and, therefore, high-energy longitu- 
dinal fiuctuations of the local spins are large and relevant. As a consequence, RPA theory 
predicts an instability of the collinear ferromagnetic ground state against some spiral order. 
Our many-body approach reveals quite a different picture. In the region of stable ferromag- 
netism we find well-defined local moments. Spin excitations in these systems are transversal 
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fluctuations of tlie local moments whose length is essentially conserved in space and time. 
Consequently, the stiffness of the variational spin-wave dispersion is always positive and 
decreases for increasing interaction strength. 

The popular "frozen-magnon" approximation for the examination of spin- wave exci- 
tations in itinerant ferromagnets can be considered as a mean-field theory in which the 
conservation of the local spins is introduced by hand. In this way, the overestimation of lon- 
gitudinal fluctuations in the Hartree-Fock/RPA approach is eliminated by construction. In 
fact, it can be shown that the frozen-magnon approximation would lead to the same results 
as our variational approach if we chose a Hartree-Fock variational wave-function instead of 
our Gutzwiller wave-function. A more detailed analysis of the frozen-magnon approach will 
be the subject of a forthcoming studylHa. 
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FIGURES 




q/A" 



FIG. 1. Variational spin-wave dispersion in x direction, E^^^{{q, 0,0)), for a generic two-band 
model with n = 0.3, J = 0.2U, and the values U/eV = 7.8, 10, 12, 13.6 correspond to 
m = 0.12, 0.20, 0.26, 0.28. The lattice constant is a = 2.5 A. 

Inset: £;™''((g, 0, 0)) and E''^''{{q/V2,q/V2,0)) for m = 0.2 and m = 0.28, respectively. The 
spin-wave dispersion is almost isotropic for strong ferromagnets. 
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FIG. 2. RPA spin-wave dispersion in x direction, £'^ ((g, 0, 0)), for a generic two-band model 
with n = 0.3 and J = 0.2[/. The value U = 4.0 eV corresponds to m < n in Hartree-Fock theory 
whereas U = 5.2 eV and U = 6.2 eV correspond to the fully polarized Hartree-Fock ground state. 
Also shown is the spin-wave dispersion for U = oo. The lattice constant is a = 2.5 A. 
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